Nonlinear Elasticity of Single Collapsed Polyelectrolytes 
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Nonlinear elastic responses of short and stiflt polyelectrolytes are investigated by dynamic sim- 
ulations on a single molecule level. When a polyelectrolyte condensate undergoes a mechanical 
unfolding, two types of force-extension curves, i.e., a force plateau and a stick-release pattern, are 
observed depending on the strength of the electrostatic interaction. We provide a physical interpre- 
tation of such force-extension behavior in terms of intramolecular structures of the condensates. We 
also describe a charge distribution of condensed counterions onto a highly stretched polyelectrolyte, 
which clarifies a formation of one-dimensional strongly correlated liquid at large Coulomb coupling 
regime where a stick-release pattern is observed. These findings may provide significant insights into 
the relationship between a molecular elasticity and a molecular mechanism of like-charge attractions 
observed in a wide range of charged biopolymer systems. 

PACS numbers: 61.20.Qg, 82.35.Rs, 87. 15. He 



O 



O 

o 



> 

o 
o 



I 

O 

o 



X 
J3 



Electrostatic interactions in aqueous media can be 
controlled by changing temperature, dielectric constant 
of the solvent, and counterion valency P, In a 

strongly charged system, correlations between counte- 
rion charge fluctuations that are neglected within the 
Poisson-Boltzmann (PB) theory give rise to various dra- 
matic phenomena which apparently contradict our un- 
derstandings based on the PB theory One of 
the most familiar example is probably the condensation 
of a single DNA molecule in the presence of multivalent 
cations ]3j ^3 j which indicates the counter-intuitive 
"like-charge attractions" A variety of problems re- 
lated with electrostatic effects is being of great interest 
in soft condensed matter physics today [ij. Ionic ef- 
fects also modifies elastic properties of PE chains drasti- 
cally d Hi, IJJ, 111 E El ■ Recent micromanipulation 
experiments have shown impressive elastic responses of 
single collapsed DNA molecules 0, |23|, which signif- 
icantly deviate from those expected from the standard 
worm-like chain (WLC) model |22i]. In spite of such ex- 
perimental evidences, theoretical and numerical studies 
of this problem have been lacked so far, mainly because 
of its nonequilibrium nature combined with long-ranged 
Coulomb interactions. 

In this Letter, we report on results of the Brownian Dy- 
namics (BD) simulations on the stretching of single PE 
condensates. The present simulation covers only a salt- 
free system with a short chain strand, but effects of coun- 
terions which are strongly coupled to PEs are taken into 
account explicitly. Force-extension {f-x) curve of a PE 
chain shows the standard WLC elasticity, a force plateau, 
and a stick-release pattern as the electrostatic coupling 
increases. This trend is consistent with that observed 
in the DNA stretching experiments 0, |23| where the 
electrostatic coupling strength is controlled by the con- 
centration of the added multivalent cations. The struc- 



ture factor of condensed couterions onto a stretched PE 
reveals a formation of one-dimensional (ID) strongly cor- 
related liquid (SCL) when a f-x curve shows a 
stick-release, while such a charge ordering is absent in the 
case of a force plateau. In order to highlight electrostatic 
effects on the macroscopic polymer elasticity, we are pri- 
marily focusing on systems with very strong Coulomb 
interactions which is beyond a typical physiological con- 
dition. Furthermore the present study deals only with 
relatively short and stiff polymers and also neglects com- 
plicated structures of biopolymers l23j , direct comparison 
with any experiment is therefore impossible at this stage. 
The primary aim in this Letter is to build bridges be- 
tween the macroscopic molecular elasticity and the the- 
oretically predicted SCL picture, which may serve as a 
first step towards microscopic understandings of available 
experiments. 

In our BD simulations, the chain comprises N har- 
monically linked beads of radius a, each of which carries 
a charge —qe. The system also includes N oppositely 
charged counterions of the same valency and radius to 
keep the overall charge neutrality. The charged par- 
ticles interact via the unscreened Coulomb interaction, 
where the solvent is replaced by a background with the 
dielectric constant e. The system is placed in a cubic 
box and periodic minimum-image boundary conditions 
are imposed to keep the monomer/counterion density 
Pmon = Pion ~ 1-5 X 10~^/a^ constaut H^Ull- Long 
time dynamics of monomers and counterions are well de- 
scribed by the position Langevin equations 
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dt 



-Vr.{/[r(t)]+f,(t), 
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where 7 is the damping constant, U the potential energy 
and fi the vectorial random force acting on the particle 
i, which satisfies the correlation property {ii{t)fj{t')) = 
2jkBTdijlS{t — t'). Because we are primarily concerned 
with dense ionic globules, hydrodynamic interactions are 
not taken into account in this study [Hl^^. The particle 
radius a, the energy scale of the Lennard-Jones potential 
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FIG. 1: The force-extension curves of a PE for three different values of the coupling parameter (a) F = 0.1 (WLC), (b) F = 40 
(force plateau), and (c) F = 110 (stick-release). The dashed line in (a) is the best fit of Eq. to the simulation data. Arrows 
indicate stretching/relaxing. 



e and r — ja^ /e are chosen as the unit of length, energy, 
and time in our simulations. The potential energy is 
composed of several kinds of contributions: U — Us + 
Ub + Uc + Ulj. The first one ensures the connectivity 
of a PE. To maintain small variances of bond lengths 
around b, i.e., the chain length around L = b{N — 1), we 
use a large value of the spring constant K = 400e/a^. 
The second one Ub accounts for the intrinsic bending 
stiffness of a chain: 
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(r,+i ~ r^) ■ (r^ ~ r^-i; 
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(2) 



The intrinsic persistent length is thus given by — 
bK/ksT. The Coulomb potential Uc acts between all 
charged particles except between bonded ones. In or- 
der to prevent a collapse of monomers and counterions, a 
truncated Lennard-Jones potential U lj acting between 



all particles is also incorporated only for jr^ 



< b. 



The coupling constant F = q^ls/a controls the rela- 
tive strength of the electrostatic interaction to the ther- 
mal energy at a distance a. (The Bjerrum length Ib = 
/ ^TTskBT K, 0.7 nm in water at room temperature.) 
In the beginning of our simulation, a randomly gen- 
erated initial configuration is allowed to collapse to a 
globule and equilibrated for 5 x 10^ time steps with 
At = 1Q~^T. Then the collapsed chain is stretched by 
pulling its one end at a constant speed, while the other 
end is fixed in a position. The stretching speed Vs is 
measured by the value /i defined as = Vs/vq, where 
vq = b/rji with tr — jb"^ N"^ / Sn^ e the Rouse relax- 
ation time for a flexible polymer with N monomers at 
ksT = e [23|. In all simulations presented below, we fix 
AT = 32, 6 = 2a, L/L° = 4, /i = 0.1 and fcsT = 0.3e, while 
r is widely changed. To complete single stretch/relax 
cycle, simulation runs for more than 10^ time steps even 
for N — 32. Simulations for larger N is desirable, but 
it is much more time consuming because the total time 
steps is proportional to Ltr oc in addition to the 
0{N'^) compulations for Coulomb interactions at every 
time step. 
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FIG. 2: The effective persistence length of a PE deduced 
from its force-extension curve as a function of F. Inset is the 
same plot on a semi-log scale. 



Figure. n(a) shows the f-x curve of the PE for F = 0.1. 
The chain remains elongated in this case, and the f-x 
curve obeys the standard WLC formula: 
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where u — x/{L) is the fractional extension |22|. Even 
in this case, the persistence length Lp differs from L9 
because of the electrostatic contribution [Tj, llJ, llM 
Ir^ IT1 | . In Fig. 12 we plot the effective Lp as a func- 
tion of F, which are determined by fitting Eq. © to the 
simulated f-x curves for 0.80 < u < 0.97. For small 
F (roughly less than 3), the electrostatic repulsion be- 
tween monomers leads to a stretched configuration of 
the chain, i.e., an increasing stiffness. As the coupling 
increases, however, it is more and more screened out by 
the counterion condensation. Far from the onset value of 
the Manning condensation (Fc ~ 2 for infinitely long, stiff 
chain H^l), the electrostatic interactions are significantly 
reduced, and Lp again decreases to Lp. Such behavior 
is also consistent with that of the radius of gyration Rg 
observed in other simulations |2^ |2^ , if one recalls that 
the persistent length is related to the radius of gyration 
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FIG. 3; Snapshots of a PE during stretching for (a) F = 40 
and for (b) F = 110. Corresponding f-x curves are Fig. Q 
(b) and (c). Dark and hght spheres represent counterions and 
monomers, respectively. For these coupling strengths, all the 
counterions are condensed onto the PE. 



as Rg ~ LpL for a semiflexible polymer. The sudden 
drop of Lri at F = 35 may be the sign of intramolecular 
collapse p^ . 

For F = 40 [Fig. ^ (b)], a distinct force plateau is 
emerged in the f-x curve, reflecting intramolecular at- 
tractions which are purely electrostatic origin. To eval- 
uate the magnitude of the condensation force between 
monomers, we calculate the excess work AW = /^^^ (/ — 

fo)dx — f{x2 — xi) from the data, where /o is the WLC 
force calculated by Eq. (PJ, and xi and X2 are chosen 
as 0.3(L) and 0.8(L}. We obtain AW - OMsTL/a, 
which roughly corresponds to / ~ LGfe^T/a. Tak- 
ing into account that / is a certain increasing function 
of F, let us tentatively compare our //F with that of 
available experiments. Assuming DNA's linear charge 
density v = q/{2a) and the charge valency g = 4 in 
our model, we have a ~ 0.34 nm (Ibp), which gives 
f/T - 0.04fcBT/bp for F = 40. On the other hand, 
for DNA with trivalent cation {q = 3) at room temper- 
ature, one finds F = 2vlB/a = 24. Using the experi- 
mental value AW = 0.33A:Br/bp for trivalent cation Co- 
Hex obtained by Baumann et al. |l9l |. one would obtain 
//F ^ O.OlfcsT/bp for F = 24, which agrees with our 
value. 

Several snapshots of the PE during stretching are dis- 
played in Fig. 121 (a) . When the tension is applied exter- 
nally, the ionic condensate (the PE-counterions complex) 
can change its overall structure to a more energetically 
stable one, because intramolecular attractions are not so 
strong (typically a few ksT as inferred above) for the 
present moderate F. Therefore, the PE can change its 
configuration continuously to minimise the stored elas- 
tic energy during stretching/relaxing, which we consider 
leads to the appearance of the force plateau. 

In Fig. n (c), we show the f-x curve for the largest 
coupling parameter: F = 110. The curve exhibits a pro- 
nounced sawtooth pattern, which may reflect intermit- 
tent plastic deformations of the glass-like ionic conden- 
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FIG. 4: Force-extension curve averaged over five independent 
runs starting from the same initial collapsed state for F = 110. 



sate [23 ■ (As will be described below, the condensed 
counterions are still in a liquid phase, but their individ- 
ual motions are extremely slow). Figure |31(b) shows a 
series of the corresponding snapshots of the PE, which 
clearly demonstrates the coil-globule coexistence within 
the single PE under the imposed extension 29] . This is a 
noticeable difference with the case in an external electric 
field where an unfolding of a PE is rather abrupt [2^ . 
For sufficiently large F, the electrostatic attractions are 
so strong that the PE condensate cannot alter its struc- 
ture easily. Then the elastic energy is more and more 
stored in the chain as the extension increases, until the 
PE condensate is locally unravelled in rather abrupt way. 

It is important to note that there is marked hysteresis 
in the f-x curve [Fig. ^ (c)] during stretching and relax- 
ing. This suggests that folding and unfolding transitions 
of the complex are first-order-like, namely, they do not 
take place at the same extension |2^. This feature is 
more distinctive in the f-x curve [Fig. 0] averaged over 
several independent runs with the same initial collapsed 
state for F = 110. This averaged curve is smoothed out 
in comparison to a single-pass sawtooth pattern as Fig.Q] 
(c) , which again indicates that each stick-release response 
is caused by a nonequilibrium configurational change of a 
PE condensate. Averaging over a large number of single- 
pass f-x curves, one would obtain a smooth f-x curve 
with vanishing peaks. However a hysteresis would re- 
main even in that limit, as long as a time scale of Vs is 
faster than the longest transition time of a PE condensate 
to overcome an intramolecular local energy barrier [s^l . 

We now turn our attention to the counterion charge 
distribution along a highly stretched PE chain. In fig. |S1 
ID structure factor S{q\\ ) is shown for F = 40 and 
110, with the PE fractional extension u kept nearly 
0.9. Here S{q\\) is calculated according to S{q\\) = 

TT Ei,j(exp[«g||(a;i(t) -a;j(i))]): where Xi{t) is the x coor- 
dinate of the i-th counterion and the bracket represents 
long-time average (taken as roughly 10^ time steps). One 
can find a crystalline order of condensed counterion for 
F = 110, while a gaslike phase is found for F = 40. The 
most prominent peak at q||a ~ tt for F = 110 is identified 
as the Coulombnic ground state configuration of the con- 
densed counterion. However, a self-diffusion of a counte- 
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FIG. 5: ID structure factor S'(gii) of counterions condensed 
onto a stretched PE with its fractional extension u kept nearly 
0.9 for r = 40 (broken hue) and for F = 110 (solid line). Inset: 
The self-diffusion of a single condensed counterion for F = 
110, which exhibits sub-diffusive behavior given by {5x^(1)) oc 
with z = 0.85 for t/r > 10\ 

rion (averaged over all particle motions) is sub-difFusive, 
though extremely slow, even for F = 110 (see the inset of 
Fig. [SI . The condensed counterions therefore form a ID 
SCL at the surface of a PE. We now expect as follows. A 
long-ranged attraction mediated by long-ranged charge 
fluctuations 11, 28] gives rise to the force plateau in f-x 
curve; on the other hand, the highly nonequilibrium elas- 
ticity such as the stick-release for large F mainly comes 
from the short-ranged strong attraction created by the 
structural correlations of ID SCL. This attraction is pre- 
dicted to dominate over the long-ranged one at a very 
short distance that is typically within the inter-particle 
distance b Isil. Because the counterions are much more 



strongly bound to the PE for F = 110 than for F = 40 
(see fig. El, it is reasonable to consider that different 
mechanisms of a net attraction between similarly charged 
monomers are underlying for F = 40 and 110. To ver- 
ify this idea, however, fine data for larger N systems are 
required because we have to diminish finite-size effects 
and influences of system-size dependent duration time of 
metastable states as much as possible. 

We expect that a sufficiently large F regime is the- 
oretically accessible from the first principle using the 
recently developed strong-coupling theory 0, 1^3 , while 
for intermediate F regime, a phenomenological approach 
which postulates a "short-ranged" intramolecular attrac- 
tion can be still useful to understand a macroscopic elas- 
ticity of a long PE [23|- (Note that the intramolecular 
attraction for moderate F is long-ranged, as mentioned, 
compared with a length scale a, but is still significantly 
"short-ranged" compared with an overall polymer con- 
tour length which often extends up to several ^m for 
biopolymers such as DNA.) 

In summary, we have studied force responses of col- 
lapsed short PE strands by the dynamic simulations in 
a wide range of F. The simulation results arc consistent 
with the experimental trends, and the idea is proposed 
by analysing the simulation data that the SCL state of 
the counterion may be responsible for the stick-release re- 
sponse. To unambigiously conclude these findings, more 
work for larger N systems will be required. 
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